Phonons in a nanoparticle mechanically coupled to a substrate 



Kelly R. Pattern and Michael R. Gcllcr 
Department of Physics and Astronomy, University of Georgia, Athens, 

(Dated: February 1, 2008) 



Georgia 30602-2451 



The discrete nature of the vibrational modes of an isolated nanometer-scale solid dramatically 
modifies its low-energy electron and phonon dynamics from that of a bulk crystal. However, 
nanocrystals are usually coupled — even if only weakly — to an environment consisting of other 
nanocrystals, a support matrix, or a solid substrate, and this environmental interaction will modify 
the vibrational properties at low frequencies. In this paper we investigate the modification of the 
vibrational modes of an insulating spherical nanoparticle caused by a weak mechanical coupling to a 
semi-infinite substrate. The phonons of the bulk substrate act as a bath of harmonic oscillators, and 
the coupling to this reservoir shifts and broadens the nanoparticle's modes. The vibrational density 
of states in the nanoparticle is obtained by solving the Dyson equation for the phonon propagator, 
and we show that environmental interaction is especially important at low frequencies. As a probe of 
the modified phonon spectrum, we consider nonradiative energy relaxation of a localized electronic 
impurity state in the nanoparticle, for which good agreement with experiment is found. 

PACS numbers: 63.22.+m, 78.67.Bf 
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I. INTRODUCTION 

There is currently great interest in properties of 
nanometer-scale mechanical systems, such as cantilevers, 
nanoparticles, and resonators .□ Because of the extremely 
small size and volume-to-surface ratio of these systems, 
interactions with their surroundings can dramatically al- 
ter their properties. In particular, it is well known that 
the vibrational spectrum of an isolated nanometer-scale 
crystal, being discrete, is qualitatively different than that 
of the same bulk material, leading to important changes 
in any property dependent on the phonon density of 
states (DOS). The differences between the vibrational 
DOS in a nanoparticle and a bulk solid are most evident 
at low frequencies: A spherical nanoparticle with diam- 
eter d and characteristic bulk sound velocity v can not 
support a mode with frequency less than about 2irv/d. 
Thus, an acoustic "gap" in the low-energy phonon spec- 
trum is present in contrast with that of the bulk, which 
has a continuous spectrum down to zero energy. How- 
ever, mechanical interaction with the environment will 
modify the discrete nature of the modes. 

In an interesting experiment by Yang et ai, El the 
phonon DOS deep inside this gap was measured in in- 
sulating Y2O3 nanoparticles. The experiment consisted 
of nanoparticles whose sizes ranged from 7 to 23 nm in 
diameter and was performed by measuring the nonra- 
diative lifetimes of an excited electronic state of a Eu 3+ 
dopant. The lowest supported mode or Lamb mode for 
nanoparticles of these sizes is about 10cm _1 . At 3 cm -1 
the DOS measured was more than 100 times smaller than 
that of bulk Y 2 3 (at 3 cm" 1 ). 

In this paper we propose and investigate a mechanism 
that could be responsible for the observed broadening 
of the nanoparticles' phonon modes. Several possible 
broadening mechanisms could be responsible for the ob- 
served effect. For example, anharmonicity leads to broad- 
ening and, therefore, to a low-energy DOS, but anhar- 



monicity is ineffective at low energy and was. found to 
be too small to account for the experiment .El Another 
possibility could be adsorbed "dirt" on the outside of 
the nanoparticle. This might lower the Q factor of the 
nanoparticle, regarding it as a resonator, reflecting a 
broadening of the vibrational modes. A third, and in our 
opinion more likely mechanism, follows from the realiza- 
tion that these nanoparticles are not isolated, but instead 
were prepared in a powered form. Thus, each nanopar- 
ticle is in weak contact with a cluster of other nanopar- 
ticles. Because the cluster is relatively large, including 
several hundred nanoparticles, each nanoparticle is me- 
chanically coupled to a reservoir that has a continuous 
vibrational DOS at low energy. This interaction broad- 
ens the modes and allows phonons in the nanoparticle to 
escape and be absorbed into the cluster. We investigate 
the affect this mechanical environmental interaction has 
on the nanoparticle's phonon DOS. 

Because we are only interested in determining the cor- 
rect origin of broadening, and do not hope to be able 
to exactly reproduce the experimental results of Ref. 
we propose the following simplified model: The cluster 
of nanoparticles is replaced by a semi-infinite elastic sub- 
strate, and one nanoparticle is placed in weak mechan- 
ical contact with it. The weak contact is imagined to 
be a few atomic bonds or small neck of material, which 
we model by a harmonic spring. For simplicity, we take 
the substrate and the nanoparticle to be made out of 
the same isotropic elastic material. Because we are in- 
terested in the low energy regime, continuum elasticity 
theory will be used to describe the nanoparticle and the 
substrate. After defining and analyzing our simplified 
model, in Sections III through [v|, we explain in Section 
[v| how the model can be adapted to address the experi- 
ment of Ref. ||, and good agreement is obtained. 

The simple model we study is related to, but differ- 
ent than, models used to stud;*, energy relaxation by 
molecules adsorbed on surfaces. □ However, in surface 
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science the interest is usually in the relaxation of rigid 
translational motion, rotational motion, or simple inter- 
nal vibrations of adsorbates. In contrast, we investigate 
the broadening of complex internal vibrational modes of 
much larger objects (which are crystalline). Also, our 
work has, much in common with that of Gurevich and 
Schober,Q where many of the same considerations and 
modeling were used to study the Lamb-mode decay rate, 
of nanoparticles caused by both anharmonicity and cou- 
pling to an enviroment of other nanoparticles. 

II. NANOPARTICLE AND SUBSTRATE 
MODEL 

As mentioned in the introduction, the model we study 
is that of a single nanoparticle in weak mechanical con- 
tact with a semi-inhnite substrate. Linear elasticity the- 
ory will be used to describe the phonons of this system. 
We assume the nanoparticle and substrate to be made 
of an isotropic non-polar material. Because we take the 
nanoparticle and substrate to be made of the same mate- 
rial, we will use the same density p and Lame coefficients 
A and fi for both. The Lagrangian for the entire system 
is given by 

L = [ d 3 r [\p(d t Mf - \\ul-pul] , (1) 

where u(r, t) is the displacement held, and 

Uij = (diUj + djUi)/2 (2) 

is the strain tensor. V is the volume of the nanopar- 
ticle, substrate, and connecting material, as shown in 
Fig. |l|. Because the Lagrangian density is local, the inte- 
gration volume in (Q) can be split into three independent 
parts: the nanoparticle, the substrate, and the connect- 
ing region. In the limit of weak coupling (diameter d c 
of connecting region much smaller than d), the surface 
area on the nanoparticle and substrate over which the 
actual boundary conditions differs from stress-free condi- 
tions are negligible, and the Hamiltonian can be written 
as 

H = -ffnano + -^sub + SH, (3) 

where H nano is the Hamiltonian for an isolated nanopar- 
ticle (with stress- free boundaries), H su b is that for an 
isolated substrate, and SH is the interaction between the 
two. The connecting region is taken to be a few atomic 
bonds or small neck of material, as shown in Fig. [I]. We 
further approximate this mechanical coupling to be a har- 
monic spring potential 

^=^:[< ano (r )-< ub (r )] 2 :, (4) 

where K is an effective spring constant, and u^ ano and 
w sub are the z components of the displacement field of 




FIG. 1: (a) Model of nanoparticle, substrate, and connecting 
region, (b) Expanded view of connecting region with dimen- 
sions d c and l c 



the nanoparticle and substrate at the point of contact, 
r . We take the z direction to be along the upward point- 
ing normal to the substrate surface. The Hamiltonians 
we will introduce below for H nano and H su ^ arc normal- 
ordered; therefore, it is necessary to normal order 5H as 
well. This operation is denoted by the colons in Eq. (Q). 

Our analysis will require the vibrational normal modes 
and spectra of the isolated nanoparticle and semi-infinite 
substrate, calculated with stress-free boundary condi- 
tions. The long-wavclcngth modes of interest here may 
be obtained from elasticity theory, to which we now turn. 



A. Isolated Nanoparticle 

Here we derive the normal modes of an isolated elastic 
sphere. The method we shall use is different than (bui 
equivalent to) that used in the classic paper by Lamb,El 
but is better suited for our purposes. The equation of 
motion of (|l|) is 

S t 2 u-«j 1 2 V(V.u)+j; t 2 VxVxu = 0, (5) 
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where v\ = y(A + 2/i)/p is the bulk longitudinal sound 
velocity and v t = \J fi/p is the transverse velocity. To 
solve Eq. (||) the displacement field can be decomposed 
into longitudinal and transverse parts, 



where 



and 



U = U[ + u t , 



V X ui = 



V • u t = 0. 



(6) 
(7) 
(8) 



With harmonic time dependence, the equation of motion 
Eq. (||) then separates into two vector Helmholtz equa- 
tions for the longitudinal and transverse parts, 



and 



fV 2 +p 2 )ui 



<? 2 )u t 



0, 



0. 



p = LU jv\ 



u>/v t . 



(9) 



(10) 



The longitudinal equation (j^) can be solved by introduc- 
ing a scalar potential 



xii = Vct) [p \ (11) 

where 4>^ is a solution of the scalar Helmholtz equation 
(V 2 + p 2 ) <fffl = 0. The transverse equation jio| ) has two 
linearly independent solutions, Ut = M and N, where 

M = V(f) (q) x r (12) 

and 

N=±VxM. (13) 

Here ^ is a solution of (V 2 + q 2 )^ = 0. The pref- 
actor 1/q is included for dimensional convenience. The 
scalar Helmholtz equations are separable in spherical co- 
ordinates and the solutions can be written as 



<^m(r) = ji(kr)Yi m (9, tp), k=p,q 



where 



(14) 



(15) 



is a spherical Bessel function of the first kind (regular at 
origin) and 



Y lm (0,<p) = (-!)' 



121 + 1 {l-m)\ 
47r (l + m)l 



P im (cos 9) e 



Here 



Plm(x) = (1 - X 2 ) 1 



dx r 



Pi{x), 



(16) 



(17) 



where Pi(x) are Legendre polynomials. 

Now we use the <pi m to construct three linearly inde- 
pendent solutions of (|J), 



Lim = ^V0; m (pr), 

M lrn = V4>i m (qr) x r, 
N ir „ e IVx M lm . 



(18) 
(19) 
(20) 



The general solution is a linear combination of L; TO , M; TO , 
and N irn , 

u(r) = ^2 [a-lrn'Llm + k m Mi m + ci m Ni m ] . (21) 



Although they are linearly independent, the vector fields 
L/ m , M/ m , and N; m are not orthogonal in space. How- 
ever, they can be written in terms of orthogonal vector 
spherical harmonics P/ m , B; m , and C; TO , defined as 



Pi m (n) = Y lm (Q)e r 



(22) 



Bi m (n) = d e Yi m (n)e e H -— — 



(23) 



C lm {Q) = -j== [—0-ee - d e Y lm (n)e v , 

(24) 

with the following properties 

^ X L ' X i'm' = Sll'Smm' (25) 

for X G B,C,P and 



dCl XL ■ X;, m , = (26) 



for X^X'. Expressed in terms of vector spherical har- 
monics, L/ m , M/ m and N; m are given by 



Urn = j[{pr)V lm {n) + ^ l< " l + 1 ^ ji(pr)B lm (n), (27) 



pr 



M lrn = y/l(l + l)Mqr)C lm (Sl), 



(28) 



and 



N,„ 



-ji(qr)Pi m {il) -+ 



qr 



qr 



ji{qr) +qrj' l (qr) B Zm (fi) 



(29) 



where prime denotes differentiation with respect to the 
argument. 

Next we impose stress-free boundary conditions 



(TijTlj = 



(30) 
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at the surface r = R of the nanoparticle. Here n is 
an outward pointing normal vector and <xy is the strain 
tensor. In an isotropic elastic continuum, 



aij = A ( V • u) Sij + 2/i ? 
In spherical coordinates ( ^p| ) implies 

a rr = 0~g r = a vr = 0. 



The three conditions (32) require that 



A(V • u) + 2/j,u rr = 0, 



U0r 



and 



In terms of the displacement field, 

uer = \{d r u e - \u e + ^d u r ), 



2 V r sin B "V 



8 v U r + d r Uu 



(31) 
(32) 

(33) 
(34) 

(35) 



(36) 
(37) 
(38) 



The boundary condition equations (32) then become 
aim. [ - Xpji(pR) Y lm + 2[ipj"(pR) Yi m ] 



and 



where 



D 
E 
F 



ci m 2fj,l(l + l)EY lm = 0, 
a im 2 Dd g Y im + h m im E Y im esc 9 

+ C\ m FdgY hn = 0, 



ai m 2im D Y lm esc 9 + h m Ed s Y lrn 
+ cimim F Y tm csc 9 = 0, 

j'i(pR) jiipR) 

R P R 2 ' 
j[(qR) j t (qR) 



R 



qR 2 



(39) 
(40) 

(41) 

(42) 
(43) 



Q3t (qR) h ^r-ji(qR) 55—- (44) 



qR 2 y qR 2 

Finally, we rewrite (|39| ) though (^lj) in matrix form as 



-Xpji(pR)Y lm + 2np3['(pR)Y lr 
2D6 6 Yi m 
2imDYi m csc 9 



2ixl(l + l)EY Zm 

imEYi, n csc 9 Fd g Y lm 
Ed 9 Y lm imFY lm csc 9 




0. 



For a nontrivial solution of ([39|) through ([IT]) to exist, the 
determinant of the above matrix must vanish. Taking the 
determinant and simplifying we find 

[ - \p :]l (pR) + 2fipj'/(pR)] E F - 4/xZ(Z + 1)E 2 D = 0. (45) 

This implies that cither 

c _fi(qR) ii(qR) 



R 



qR 2 







(46) 



or 



Xpjt ( P R) + 2 m '/(pR)] F — 4pl(l + 1)E D = 0. (47) 
If (|4^) is met, then this imposes certain constraints on 



aim Mm and q m , which require ai r 



0. This can 



easily be seen in the above matrix by setting E = 0. If 
( [47| ) is met, bi m has to be zero. In conclusion, we have 
two branches of vibrational modes. The branch in which 
([4(f) is satisfied, 



u(r) 



>(r), 



(48) 



are referred to as the torsional modes, where n specifics 
the radial quantum number. The other branch is found 



I 

when ( [47|) is satisfied, 

u(r) = ai mn Li mn (v) + c; m „N im „(r), (49) 

which are the spheroidal modes. 

To quantize the. vibrational modes we write the dis- 
placement field asQ 



Unano(l") ^ 



2pU!,J 



Qj *Mr)+4*}(r) , (50) 



where 



J = [S or T, n, 1, m] 



(51) 



is a label uniquely specifying a nanoparticle cigenmodc. 
The first entry S or T specifies whether the mode is in the 
spheroidal or torsional branch, respectively n is the ra- 
dial quantum number and I and m are the usual angular 
momentum quantum numbers, a and are phonon an- 
nihilation and creation operators which satisfy the Bose 
commutation relation 



[a j, a?j,} 



TJ.J> 



(52) 
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The ^ j arc vibrational eigenvectors normalized such 
that 



/ \*j(v)\ 2 d 3 r = l, 
Jv 



(53) 



where V is the volume of the nanoparticle. Assuming 
(without proof) that the modes \& j form a complete set, 



(54) 



it can easily be shown that u satisfies the correct equal- 
time canonical commutation relation with tt = pdtu, 
namely 



[u l (r),7rV)] =ift* tf *(r-r / ). 



B. Isolated Substrate 



(55) 



The vibrational modes for a semi-infinite isotropic elas 
tic substrate, with a free surface at the xy plane and ex 
tending to infinity in the negative z direction, were quan 
tized previously by Ezawap therefore, the details will 
left out here. The displacement field can be written ai 



u S ub(r) = 



2pujj 



bMr) + b]ff(r) 



(56) 



where b and b^ are the annihilation and creation oper- 
ators for the substrate phonons. The index I, like the 
index J for the nanoparticle, uniquely specifics a phonon 
mode for the substrate, fj are eigenfunctions of ([|) sub- 
ject to stress-free boundary conditions at the z = plane. 

In what follows we will need the spectral density of the 
isolated substrate which is defined as 



iV su b(r,w) 



1 



Im^subOMNw) 



(57) 



where D^.ir, r', u>) is the Fourier transform of the re- 
tarded phonon Green's function 



sub V 



(r, r', t) = -i6(t) ( [< ub (r, t), < ub (r', 0)] ) (58) 

of the substrate. The spectral density at the free surface 
of silicon, regarding it as an isotropic elastic continuum, 
was calculated in Appendix B of Rcf. [| There we ob- 
tained 



iV S ub(w) = CsiW, C Si 



1.4 x 10~ 46 cm 2 s 2 . (59) 



III. GOLDEN-RULE LIFETIMES 

The relaxation rate or inverse lifetime of the perturbed 
eigenmodes of the nanoparticle can be calculated using 
Fermi's golden rule (setting h = 1), 



2 7 r]T|(f|^|i)| 2 ^ i -c, f ), 



(60) 



where the initial and final states are 

|i)=4l°) and |f)=&/|0>- 
Using Eqs. (@), ©, and © leads to! 



7r^ 2 |*}(r„)| 2 _|//(r ) 



J 2p 2 
Noting that 



E 



-6(wj 



(61) 



(62) 



J2 \f?(ro)\ 2 5(u.j - wj) = 2 P LjN snh (cj), (63) 



we obtain (reinstating factors of K) 



irK 2 N suh (ujj) 

hp 



*}(ro)| 5 



(64) 



Using K = 1.0 x 10 4 erg cm -2 and p = 2.3 gem -3 , which 
are appropriate (see Ref. ||) for a weak link in Si, relax- 
ation rates and Q factors are given in Table [j] for some low 
lying modes. The Q factor is defined here as r divided 
by the period T, 



^ T 2tt 7 ' 



(65) 



where 7 = Tit -1 is an energy width. 

The values of the Q factors we obtain for the low lying 
modes are incredibly large, reflecting the fact that the 
reservoir (substrate) is extremely ineffective at absorb- 
ing energy at these low frequencies. As we will discuss 
below in Section VI, the lifetimes (and Q factors) for the 
model considered here cannot directly be compared to 
the experiment of Rcf. ^| without accounting for the dif- 
ference in sound speeds between a solid Si substrate and a 
weakly bound nanoparticle cluster, as well as some other 
less important modifications. There we shall show that 
the coefficient C in Eq. (|5^) should be enhanced by about 
1000 before making such a comparison, which decreases 
the Q factors by this same factor. However, the Q factors 
corrected in this way are still huge, and the good agree- 
ment with the observed low-frequency DOS (see below) 
suggests that the Q factors of the nanoparticles studied 
experimentally in Ref. ^ are also very large. 

In Ref. |l0| we used the golden rule result (Q) to esti- 
mate the phonon DOS at low energies. This is achieved 
by replacing, in accordance with Fermi's golden rule, 
each discrete mode in the isolated nanoparticle by a 
Lorentzian with a width given by d6~i|). (More precisely, 
this amounts to approximating the energy-dependent 
phonon self-energy for each mode J with its value at 
u) = u>j, a procedure often called the quasiparticle-pole 
approximation.) However, this procedure is unreliable at 
low energies because the actual line-shapes of the broad- 
ened modes are non-Lorentzian in the tails. Nevertheless, 
we obtained a DOS at 3 cm. -1 that was only 20 times 
smaller than that observed .Ell 
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TABLE I: A few representative relaxation rates and Q factors. 
Note that only the m = spheroidal modes are broadened by 
the interaction. 



(S,Z,m,n) 


u) (rad s 1 ) 


r-V 1 ) 


Q factor 


(S,2,0,l) 


3.5 x 10 12 


2.5 x 10~ 4 


2.2 x 10 15 


(S, 1,0,1) 


4.9 x 10 12 


1.0 x 10 -8 


7.8 x 10 19 


(S, 0,0,1) 


1.1 x 10 13 


3.6 x 10" 4 


4.9 x 10 15 


(S,0,0,2) 


1.7 x 10 13 


1.1 x 10" 3 


2.5 x 10 15 



IV. MANY-BODY THEORY OF THE DOS 
A. Local DOS 

To leading order in the elcctron-phonon interaction 
strength, the electronic population relaxation rate due 
to phonon emission (for example, as measured in Ref. |^) 
is given by Fermi's golden rule, which states the rate (for 
a deformation potential electron-phonon interaction) is 
proportional to the square of the electron-phonon cou- 
pling strength times the phonon DOS. In a translation- 
ally invariant system the DOS does not have any position 
dependence, but in a nanoparticle one must distinguish 
between the "global" DOS (the DOS relevant for ther- 
modynamics) and the local eigenfunction-weightcd DOS, 
which is the one that determines phonon emission rate. 
We will call this position-dependent DOS the local DOS, 
and denote it by g(r, ui). The precise definition of g(r,u>) 
will be given below. 

From a theoretical point-of-view, the quantity describ- 
ing the local vibrational dynamics in the nanoparticle is 
the (retarded) phonon Green's function 

D%(r,r',t) = -ie(t)([u i (r,t),u j (r',0)}) H , (66) 
where 



_ Tr(e-^Q 
[ >H ~ Tte-f* H ' 
and with the Hamiltonian given by 

H = H Q + SH. 



(67) 



(68) 



Here Hq is the Hamiltonian of the isolated nanoparticle 
and substrate, 



S ' LUja/jdj 
j 



i 



(69) 



and, as mentioned in Section II, 5H is a harmonic spring 
potential given in Eq. ([!]). 

In this section the phonon Green's function D always 
refers to the nanoparticle, and the label "nano" will be 
suppressed. The imaginary part of the Fourier transform 
of Dp (r, r, t) defines the nanoparticle's phonon spectral 
density 



1 



Im fl^(r,r,w). 



(70) 



For an electron system (or any system of particles), the 
spectral density defined above is precisely the local DOS. 
However, because the elasticity equation of motion (||) 
is second order in time, the spectral density and DOS 
(both local and global) differ by a factor of 2pu. In ad- 
dition, the vibrational spectral density (|7^) is a tensor, 
whereas the phonon emission rate probes some coupling- 
constant-weighted sum of tensor elements. Because we 
are ascribing the observed reduction in phonon emission 
(in going from bulk to nanoparticle) to a reduction in 
the local DOS, our results are not sensitive to the pre- 
cise way in which a scalar quantity is constructed from 
the tensor, as long as the same measure is used in both 
the nanoparticle and bulk. It will be most convenient to 
investigate the trace of the local DOS tensor. Therefore, 
the quantity we calculate in this paper is 



2pw N * 



(71) 



which we shall refer to as the local DOS. g(r,ui) char- 
acterizes the number of states per unit energy per unit 
volume near position r. In a bulk material with Debye 
spectrum, (f7lf) reduces at low frequency to 



f(r,w) 



2^ 



1 



(72) 



independent of r. Eq. ( |72[) is the well-known Debye for- 
mula for the vibrational DOS of a crystal. 

The local DOS g(r,u>) controls the phonon emission 
rate for an impurity atom sitting at position r. Although 
the impurity locations in a real nanoparticle are assumed 
to be random, dopants near the surface are known to 
be optically inactive; hence, the experiments (including 
that of Ref. |^) do not probe the phonon DOS near the 
nanoparticle surface. Therefore, we introduce a particu- 
lar volume-averaged DOS 



g{b,w) = 



J d 3 rg(r,u;) 

r<b 



H 3 



(73) 



which characterizes the average <j(r, w) within a sphere 
of radius b. In the limit b — > R, in which case the local 
DOS is averaged over the full nanoparticle volume, we 
obtain the global (or thermodynamic) DOS, which, for 
an isolated nanoparticle, is 



g(R,cj) = -^5(, 



(74) 



Physically, we expect b to be somewhere between R/2 
and R. 



B. Perturbative Calculation of the Local DOS 

The retarded Green's function ([36j) for the nanoparti- 
cle can be obtained by calculating the Euclidean time- 
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ordered (or imaginary time) Green's function defined by 
/>>,r',r) = - (7V(r, r)u^r', 0)) H . (75) 
In the interaction representation, 



DOS of Nanoparticle 



Tu 



{v,T)ui{v\Q)e-^ SH{T ' )dT '\ 

I Ha 



- f SH<T')dT 
I Jo v ' 



H 



(76) 

where the expectation values are with respect to Hq. By 
expanding the exponentials to leading order in the per- 
turbation and Fourier transforming, ( |7q ) can be written 
as 



I>«(r,r» = ^(r,r>) + E / D o^ r " 

7,/ ** 



x n w (r", r'", w)D$ (r'", r', w) dV'd : 



A-/ 



(77) 



where 



1 



1 



D « (r , r >) = E 2M£)fM 

^y> Zputj [ioj — uij ioj + ujj 

(78) 

is the free propagator and IT-' is the, leading order self- 
energy, given at zero temperature byu 



2 r 



1 



1 



x 5 iz 5 jz S(T-r )S(r' -r ). 



(79) 



The fi(r) are the substrate eigenfunctions discussed in 
Section II, and ro is the point at which the nanoparticle 
is connected to the substrate. Retarded quantities are 
obtained by analytically continuing iu — > u) + i0 + . 

To calculate the local DOS ( fn] ) we need to solve the 
Dyson equation for the nanoparticle Green's function, 
written symbolically as 



D = Dn + DnUD. 



(80) 



A solution to (|30|) can be obtained by introducing matrix 
representations for D, Do, and II, in which case 

D = (D„ 1 -IT)- 1 . (81) 

The matrix representation we use is defined by 

0(J,J',lu) = Y^ J d 3 rd 3 r'¥j(r)y j j,(r')O ij (r,r',u), 

ij V 

_ (82) 
where O = D, Dq, or II. In ( p2[ ) the integration is over 
the volume V of the nanoparticle, and the *ffj(r) are the 
nanoparticle eigenfunctions. The inverse transformation 
is 



0« (r, r', w) = ^ ¥j(r) 0(J, f, w) ^(r') 



(83) 



n 0.6 - 
3 



O 
D 



Energy (wavenumber) 

FIG. 2: The phonon DOS, given in states per wavenumber 
per cm 3 , of a 10 nm Si nanoparticle, weakly coupled to a 
semi-infinite substrate. 



A nanoparticle with a diameter of 10 nm has approxi- 
mately 8,000 atoms in it, thus, there are roughly 24,000 
acoustic vibrational modes. By knowing the number of 
modes, a Debye energy can be defined: the Debye energy 
is the energy at which there are 24,000 elasticity-theory 
modes that lie below in energy. For our nanoparticle, the 
Debye energy is about 240 cm -1 . The Debye energy cut- 
off truncates the Hilbcrt space, which leads to finite-size 
matrices.Ej This enables every mode J of the nanoparti- 
cle to be included in the calculation of the Green's func- 
tion 



V. RESULTS 

In this section we present our results for the phonon 
DOS in a 10 nm Si nanoparticle, obtained by solving 
the Dyson equation ( pO] ) for the phonon Green's func- 
tion, as explained above. As we have discussed, the DOS 
g(r,uj), defined in (f7l|), is a local quantity that varies 
with position within the nanoparticle, and, as mentioned 
in the previous section, the quantity we are interested in 
is g(b,u>), which is <?(r, to) averaged over a sphere of ra- 
dius b centered at the center of the nanoparticle. Because 
we have found no significant dependence of g(b,u>) on b, 
for the physically relevant values of b (R /2 < b < R), we 
plot g(R,uj). As stated above in Section |1V A , g{R,u>) is 
the global phonon DOS in the nanoparticle. 

For simplicity we assume both the nanoparticle and 
the substrate to be made of Si; this allows us to use the 
surface spectral density (^) calculated in Appendix B of 
Rcf. ||, where Si is treated as an isotropic elastic contin- 
uum with longitudinal and transverse sound velocities 



v\ = 8.5 x 10 5 cms \ 



././' 



5.9 x 10 5 cms -1 , 



(84) 
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FIG. 3: Vibrational DOS at low energies. 
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FIG. 4: Expanded view of the low-energy DOS. Note that 
the DOS vanishes at zero energy, as expected. 



and mass density p = 2.3 g cm -3 . In the final section of 
this paper, where we compare our results to the experi- 
ment of Ref. ^ we will introduce an important correction 
to account for the differences between a solid Si substrate 
and a nanoparticle cluster. 

In Fig. H the global DOS g(R, uj) of a 10 nm diameter 
nanoparticle is given up to 100 cm" 1 . The modes above 
100 cm" 1 were included in the calculation, but the long- 
wavelength approximation of elasticity theory becomes 
invalid at high energy. Thus, only the lower part of the 
spectrum is shown. Fig. shows the low-energy phonon 
DOS up to about 15 cm The large peak on the right 
side is the well-known Lamb mode. The phonon DOS at 3 
cm" 1 is approximately 4.9 x 10 10 states per wavenumber 
per cm 3 . 



VI. COMPARISON WITH EXPERIMENT 

In this section we compare our results with the experi- 
ment of Ref. |^, where the one-phonon emission rate (and, 
therefore, the phonon DOS at 3 cm -1 ) in a cluster of 
Y2O3 nanoparticles was observed to be 8.2 x 10~ 3 times 
that in bulk Y2O3. In particular, the excited 5 Di(II) 
state of Eu 3+ in the nanoparticles had a phonon-emission 
lifetime of 27 ps, compared with a bulk value of 221 ns. In 
order to make a comparison of our results to that of the 
experiment, two modifications of our calculation have to 
be performed. 

In our model, the cluster of nanoparticles has been 
replaced by a solid substrate. However, the spectral den- 
sity (|5^) of the substrate, which at long wavelengths, 
is determined by the sound speeds and mass density of 
Si, is much smaller than that of the nanoparticle clus- 
ter. Treating the long-wavelength modes of the cluster 
with elasticity theory (or, even simpler, approximating 
the random cluster by an ordered cubic lattice), shows 
that the spectral density (J5£j) should be replaced by (the 
subscript "cl" referring to cluster) 



whcrcl 



iV c i(w) = C C \LU, 



<-"cl — —5 <-^Si- 

V c\ Pel 



(85) 



(86) 



Here v c \ is a characteristic sound speed in the cluster 
and p c \ is its mass density. The 1 / w 3 j dependence in ( |8"6| ) 
comes from the well-known velocity dependence of the 
Debye DOS, and the l/p c \ factor comes from the defini- 
tion of spectral density [see discussion following Eq. 
Approximating the cluster by an ordered cubic array with 
lattice constant d (the nanoparticle diameter), yields 



and 



Pel ~ -^psi 
6 



(87) 



(88) 



where K is the effective sprin g co nstant connecting the 
nanoparticles, given after Eq. (p4), and 



4 ( d 



(89) 



is the mass of one nanoparticle. Using d = 13 nm (the 
mean diameter in Ref. ^), we obtain an enhancement 
factor 



Csi 



9.8 x 10' 



(90) 



This factor increases the nanoparticle DOS (at frequen- 
cies below the Lamb mode) by nearly three orders of 
magnitude. 
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There are several other marginally important correc- 
tions, most of which will be ignored, and one that will 
be included just for completeness: In the model analyzed 
above, the nanoparticle was connected to its surround- 
ings by only a single contact point, whereas a nanoparti- 
cle in a cluster most likely has more than one connection. 
As the number of contacts increases, this simply scales 
the DOS (away from the peaks) linearly with the num- 
ber of contact points. Conservatively, we expect that the 
multiple contact points present in the real system will 
increase the decay rate of the nanoparticles' modes, and 
hence the phonon DOS well below the Lamb mode, by a 
factor of two. 

The following additional corrections have also been 
considered and were found not to be significant (and are 
not included in our final results): 

(i) The actual experiment of Ref. |^ was done on an en- 
semble of nanoparticles with mean diameter 13 nm and 
standard deviation of 5 nm. To understand the effects 
of this size distribution, we have calculated the DOS at 
3 cm -1 averaged over a Gaussian distribution of diame- 
ters centered at 10 nm. Even for very wide distributions 
(standard deviation up to 8 nm), the ensemble averaged 
DOS-at 3 cm" 1 is increased by no more than a factor of 
twoJlil In addition, the correction for re-centering the size 
distribution from 10 nm to the experimentally observed 
13 nm average size, leads to corrections of order unity. 

(ii) Our calculations were done for a nanoparticle and 
substrate made of Si, while the experiment was done 
on Y2O3 nanoparticles, which, of course, has different 
mass density and sound velocities. The differences in 
sound speeds and mass only shift around the modes of 
the nanoparticle by a small amount; this has the same 
effect as a small change in the diameter of the nanopar- 
ticles, which we have found to be negligible. As for the 
substrate (or more precisely the replacement of the sub- 
strate with a cubic lattice of nanoparticles), the change 
in mass density does effect the spectral functio n (p9 |) by 
changing the velocity ( |S7l) and mass density (|88|), but 
this change is only of order unity. 

(iii) The nanoparticles of Ref. || were immersed in He, 
either liquid (for T<4.21 K) or gas (T>4.21 K). How- 
ever, the results were found not to change through the 
liquid-gas transition, presumably because of the large 
sound-speed mismatch between supcrfiuid He and Y2O3. 
Therefore, we have ignored the presence of He in our 
theory. 

(iv) The experiment of Ref. ^| was done at tempera- 
tures between 1.5 and 10 K (excluding the interval 2.17 
to 4.21 K), whereas our calculations assume zero tem- 
perature. The effect of finite temperature is to stimulate 
phonon emission into the bath (i.e. substrate or cluster). 
However, this is not important until the Bose distribu- 
tion function of the bath at the Lamb mode frequency 
(approximately 13 cm" 1 ) becomes of order unity, which 
does not occur until the temperature exceeds about 18 
K. 

Taking into account the first two modifications, and 



ignoring the rest, we obtain a 3 cm 1 DOS given by 

states 



g(R,3cm- L ) = 4.9 x 10 1 
x 980 x 2 



wavenumbcr cm 



9.6 x 10" _. (9i) 



wavenumbcr cm 



The DOS of bulk Si is given by 



E 



= 3.9 x 10^E Z 



states 
wavenumber cm 3 



(92) 



where E is the energy in wavenumbers. Thus, the theo- 
retical ratio of the DOS of the nanoparticle to that of the 
bulk material at 3 cm" 1 is 2.7 x 10" 3 . As stated above, 
the experimental ratio of nanoparticle to bulk DOS was 
found to be approximately 8.2 x 10" 3 . The agreement 
between our theory and the experiment of Ref. ^ is ex- 
cellent considering the simplicity and robustness of our 
model. We conclude that the low-energy phonon DOS 
observed in Ref. || is consistent with our cnviromental 
broadening mechanism. 



VII. CONCLUSIONS 

Motivated by an interesting experiment! measuring 
the low-energy phonon DOS in a insulating nanoparticle, 
we have thoroughly investigated a simplified model of a 
single nanoparticle weakly coupled to its environment, a 
semi-infinite substrate. The environmental interactions 
were found to significantly affect the DOS at energies 
below the Lamb mode. 

Additionally, we have used the results of our model to 
predict the effect of environmental interaction in a clus- 
ter of nanoparticles like that studied in Ref. 0. Although 
it is necessary to estimate the value of several quantities 
appearing in the model, we believe that we can do this 
accurately enough to obtain a final result that is correct 
at the order-of-magnitude level, with no free parameters. 
Because our results for the 3 cm" 1 DOS is only about 
three times smaller than that observed in Ref. our 
broadening mechanism and the resulting phonon spec- 
trum is clearly consistent with that experiment. 
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We note, however, that when the Lamb mode frequency 
becomes resonant with the electronic two-level system, 
the electronic relaxation rate is not simply related to the 
phonon DOS. 



